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03 Abstract. 

The macroscopic dynamics of a rotating superfluid deviates from that of a simple perfect 
fluid due to the effect of vorticity quantisation, which gives rise to a substructure of cosmic 
string type line defects that results in a local anisotropy whereby the effective average pressure 
in the direction of the vortex lines is reduced below its value in lateral directions. Whereas 
previous descriptions of this effect have been restricted to a non-relativistic framework that is 
adequate for the treatment of liquid helium in a laboratory context, the present work provides a 
fully relativistic description of the kind required for application to rotating neutron star models. 
To start with, the general category of vortex fibration models needed for this purpose is set up 
on the basis of a Kalb-Ramond type variational principle. The appropriate specification of the 
particular model to be chosen within this category will ultimately be governed by the conclusions 
of microscopic investigations that have not yet been completed, but the results available so 
far suggest that a uniquely simple kind of model with an elegant dilatonic formulation should 
be tentatively adopted as a provisional choice so long as there is no indication that a more 
complicated alternative is needed. 
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1 Introduction. 



Until now, calculations of the macroscopic effect of quantised vortices in the superfluid layers 
of neutron stars have relied on the use of non relativistic formalism whose most complete 
and detailed development is due to Lindblom and Mendel' 1 '. Despite the fact that relativistic 
corrections are by no means small, their neglect has so far been physically justifiable in the 
analysis of the most important observational effect for which they are relevant, namely pulsar 
frequency glitches, due to the many major uncertainties that are involved in the vortex pinning 
model that has been proposed as an explanation' 2 ' . Although there is not yet any urgent 
need for the higher accuracy that it can provide, the use of a fully relativistic formulation is 
already wanted in order to satisfy a more pressing requirement, namely that of compatibility 
with the fully relativistic fluid models which have long been in standard use for representing 
the the gross structure of neutron stars: the main drawback of the mongrel combinations of 
relativistic and non relativistic elements grafted together in the treatments that have been used 
so far is not so much their lack of physical accuracy, but rather their lack of mathematical 
coherence which makes them very unweildy to work with. The purpose of the present work is 
to overcome this disadvantage by developing a more conveniently coherent approach in which 
the macroscopic effect of vortex quantisation is allowed for in a fully relativistic framework. 
The use of this approach is exemplified here by the introduction of a simple prototype model 
that is the natural relativistic generalisation of the simplest special case within the category 
of non-relativistic models considered by Lindblom and Mendel' 1 ', allowing only for the bulk 
motion of the neutron matter and ignoring the independent motion of the residual protons, 
whose analysis is postponed for future work of a more detailed nature. 

The concepts on which the present approach is based were originally developed not in neu- 
tron star theory but in a more exotic cosmological context, specifically that of axion field theory. 
It was Lund and Regge' 3 ' who first drew attention to the analogy between Kalb-Ramond cou- 
pled string motion in an axion field and vortex motion in an ordinary non relativistic superfluid. 
The immediate inspiration for the coherent relativistic theory presented here was provided by 
the more recent analysis of Davis and Shellard' 4 " 5 ' , and its modification to allow for more gen- 
eral kinds of field by Ben Ya'acov' 6 " 7 ' . The outcome effectively extends and completes the 
results of an earlier pionneering investigation using a different but ultimately related approach 
by Lebedev and Khalatnikov' 8 ' . 



2 



2 The perfect fluid limit. 



The proposed theory is obtained as the natural generalisation of an appropriately reformulated 
version of ordinary "barytropic" perfect fluid theory given ^ by a Lagrangian variation principle 
which will now be described. 

By definition, perfect (relativistic) fluid medium is characterised by a stress momentum 
energy density tensor T^ u that is expressible in terms of its timelike unit eigenvector u^ 1 , by 
the familiar formula T^ v = pu^u u + P(g^ u + u^u"), with u^u^ = — 1 where g^ u denotes the 
background metric, and the scalars p and P are respectively the relativistic mass-energy den- 
sity and the pressure. The particularly simple case of an ordinary "barytropic" perfect fluid 
model is specifiable by an "equation of state" that determines p as a function of the number 
density n of a conserved current, or equivalently, in conjugate form, that determines P as a 
function of the relativistic chemical potential or "effective mass" per particle variable p. The 
mutually conjugate functions for the mass density and the pressure are not independent but 
are symmetrically related by a Legendre type transformation given by 

dp dP 

p{n\ + P{p\ = np , ^ = ^' n = ~dp^' ( ' 

These functions determine the corresponding sound speed, c s say, by a formula of the familiar 
form 

c 2 = = l^t. (2 2) 

s dp p dn 

In such a model, as also on the generalisation to be presented in the next section, a par- 
ticularly important role is played by the vorticity 2-form w^ v , which is defined as the exterior 
derivative 

«V = 2V^„], (2.3) 

of the relevant 4- momentum covector, p\, which is the dynamical conjugate of the conserved 
particle current vector n x . The momentum covector is related to its dual, the corresponding 
momentum trivector, by 

^ = \e^ vp E^ , (2.4) 

where S\p, vp is the alternating tensor of the 4-dimensional spacetime background, while similarly 
the conserved current vector is given in terms of its dual 3-form N^ vp by 

n X = ^ p AU . (2.5) 

The condition that the current is conserved, i.e. V M n M = 0, is expressible (in terms of covariant 
derivation with respect to the (flat or curved) spacetime metric g^ u ) as the closure of the 
corresponding 3-form, i.e. 

V [M iV, H = (2.6) 
(where the square brackets denote index antisymmetrisation). 
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In the particular case of the barytropic perfect fluid model, though not in the generalisation 
to be given in the following section, the momentum covector and its dynamical conjugate the 
current vector will be expressible in terms of the quantities introduced in (2.1) simply by 

Hx = fiu x , n x = nu x , (2.7) 

which implies 

n 

In this case the complete set of equations of motion is given by supplementing the kinematic 
relation (2.6) by a dynamical momentum transport equation which is obtainable from a con- 
vective variational principle' 9 ^ in the form 



H^ P V X H X ^ = , (2.9) 



whose more familar dual form is 



n p w^ = . (2.10) 

The most familiar category of equations of state to which this theory can be applied is 
the "polytropic" power law kind, p oc n 1 , for a constant "polytropic index" 7. Particularly 
important examples are the pressure free ("dust") case, 7=1 (for which the flow lines are 
geodesies), and the ultrarelativistic (massless particle) gas case, 7 = 4/3. In the generic case 
the variable fi will not vary proportionally to n so that the momentum 3- vector H^ vp will not be 
simply identifiable with the index raised version N^ vp of the current 3-form N pvp . However such 
an identification will be possible in the special case of the "stiff" (linear) Zel'dovich equation 
of state 1101 , with 7 = 2, as given by 

k 3 h 3 

N N^ up H — — N (2 11) 



' 12m 2 ^ ' m 

where m is a fixed mass scale (exact identification being obtainable by choosing units such that 
m 2 = ft 3 ). 

For an arbitrary equation of state there is a corresponding (zero temperature) superfluid 
model^ 11 ' included as a particular application of the foregoing formalism when the (automatically 
conserved) vorticity is zero, since by (2.3) and (2.4) we shall then have 

«V = V p H p ^ = 0, (2.12) 

which is the necessary and locally sufficient integrability condition for the existence of a gauge 
scalar ip that satisfies 

//„ = V„(M H pi/p = he>" pa V a cp , (2.13) 

and that represents a condensate order parameter which may be presumed to be axionic in the 
sense of being periodic, so that it may have string-like topological defects. The normalisation 
used here is such that periodicity 27r is consistent with the usual quantisation condition that 
the conserved circulation integral, 

k = j n v dx v , (2.14) 
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round a defect should be a multiple of the standard Bohr unit 2nh. In the special case of a 
Zel'dovich fluid the corresponding superfluid model reduces exactly to the standard massless 
(Goldstone boson) axion field model' 12 " 13 " 14 ] with the coefficient m representing the correspond- 
ing Higgs mass scale. 

3 The vortex dynamical generalisation. 

A Kalb-Ramond field is definable, modulo addition of the exterior derivative of an arbitrary 
1-form a^, for any (perfect or other) conserved fluid as a gauge 2-form Bp V from which the 
physical current 3-form is obtainable by exterior differentiation, i.e. 

N^p = 3V[ M -B„ P ] , B^y = Bp V + 2V[ M a^] , (3.1) 

the closure property (2.6) being the Poincare integrability condition that is necessary and locally 
sufficient for its existence. One can use it for setting up an appropriate action principle in which 
B^ v is to be treated as a free field with the 2-form Wp U as its source, which means that, instead 
of remaining a secondary construct obtained via (2.3), the vorticity must be promoted to the 
status of an independent - albeit not entirely free - field in its own right. This makes it eligible 
to contribute directly to the variational action in such a way as to allow at a macroscopic level for 
the effect of a smooth distribution of microscopic quantised vortices in a realistic (compressible) 
model for the treatment of superfluid dynamics on a large scale, in contexts such as that of 
a neutron star interior, as well as ordinary laboratory applications involving liquid Tie, in the 
zero temperature limit. (To deal with the Landau type 2-constituent theory that is needed 
to allow for the independent entropy flux in a superfluid at finite temperature' 8 " 11 " 15 ', a more 
intricate theory would be required.) 

The required form of the Lagrangian is given' 9 ' by 

C = A- \e^B^w pc , (3.2) 

where the "master function" , A, is a gauge independent scalar function of the current 3-form 
N^ vp (and also, of course of the background space time metric), which is specified simply by 
taking A = —p. The innovation to be considered here consists in allowing the master function 
A to depend not just on Np Vp but also on the vorticity 2-form Wp V . Consistently with the 
previous definition' 9 ' in the perfect fluid case, the current momentum trivector H IMUp and also, 
in the generic case, the relevant vorticity momentum bivector , are specified as partial 
derivatives by the (fixed background) variation rule 

5A = -±H^8Np U p - l -V v 5wp V . (3.3) 

The final coupling term in (3.2) is contrived so that requiring invariance of the corresponding 
action integral with respect to free variations of the gauge field B^, for a given value of Wp U , 
is equivalent to imposing a field equation 

VpH^ = \e up ^w^ . (3.4) 
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It can be seen that this is precisely what is needed for the on shell value of w^ v to be consistent 
with the original definition by the equations (2.3) and (2.4). 

It remains to promulgate the rule governing the specification in (3.2) of the source field w^ u , 
which obviously can not be entirely free since otherwise the gauge invariance (3.1) would be 
violated and the system would in consequence be overdetermined. To get the desired result, the 
vorticity is introduced^ via a fibration in a manner analogous to that by which the current was 
introduced in the convective variational formulation^ , but using a base space of only 2 (in- 
stead of 3) dimensions with a vorticity measure 2-form having components Wij with respect to 
local coordinates \ l sa Y f° r i = 1? 2, which act as a pair of freely variable (diffeomorphism gauge 
dependent) scalar fields. A projection i— > \ l structures the relevant spacetime neighbour- 
hood as a fibre bundle over the material base manifold with 2 dimensional fibres interpretable 
as vorticity world sheets, and the base measure form then determines the corresponding 
vorticity 2-form as its spacetime pullback: 

«V = w ijX l ^X J ,u • ( 3 - 5 ) 

It will always be possible locally to choose the coordinates to be flat with respect to the 
symplectic structure Wij on the base space, so that its non zero off diagonal components have 
unit magnitude, which gives the explicit expression 

«V = 2 X>X%] • ( 3 - 6 ) 
The prescription (3.5) automatically ensures that w^ v is closed, 

V^uvp] = , (3.7) 

which is needed both as a condition for gauge invariance of the action integral and also as an 
integrability condition for the field equation (3.4). It also automatically ensures that w^ v is 
algebraically degenerate, i.e. its determinant \w\ vanishes, so that it satisfies 

e^ pa w^w p<7 = , «v = ^we^E*" , (3.8) 

where w is the scalar vorticity magnitude and E^ v is the unit bi-vector (as normalised by 
E^ v E Vii = 2) tangential to the congruence of vortex flux 2-surfaces characterised by constant 
values of the pair of comoving coordinates x 1 - The allowed variations of are those generated 
by infinitesimal displacements of the form 5\ l = — ^dvOCi f° r an arbitrary spacetime vector 
field so they are given by the corresponding Lie derivation formula which - in view of (3.7) 
- takes the form 

<S«V = 2V [M K ]p e) . (3.9) 

Imposing that the action integral be invariant with respect to such variations leads finally to a 
basic dynamical equation given by 

uv(ra" - VpA"") = , (3.10) 
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which is in qualitative agreement with a form previously proposed by Lebedev and Khalatnikov® 
on the basis of a different approach whose relationship to the present analysis will be described 
in the appendix. 

The covariance of the foregoing formulation implies that the system will satisfy a corre- 
sponding Noether identity of the standard form 

V u T pu = (3.11) 

with 

Working out the extra term that will be induced on the right of (3.3) by a virtual variation 
5g^ v of the background metric, the relevant stress energy momentum density tensor is found to 
be given by 

Tt = \H^N vpc + \»°w Vff + kg* . (3.13) 



4 Conjugate reformulation. 



It is useful for many purposes, and in particular for relating the present approach to the previous 
work of Lebedev and Khalatnikovt 8 ' , to reformulate the theory that has just been presented in 
a dynamically conjugate version in terms of a modified master function ^ that is obtained from 
the original master function A by a Legendre type transformation that is derivable as follows. 
We start by replacing the Kalb Ramond form by its dual 

b^= l -e^B pa , (4.1) 

so that by (3.1) the current vector (2.5) is expressible simply by 

n p = y u b pv , (4.2) 
which entails the current conservation law 

Vy = , (4.3) 

as an immediate consequence. 

Instead of postulating (3.1) or its dual version (4.2) as an axiom, we can impose it via the 
variational principle by introducing an appropriate Lagrange multiplier 7r p and replacing the 
original Lagrangian C as given by (3.2) by the correspondingly augmented Lagrangian, whose 
form, after subtraction of a divergence which does not affect the field equations, is given by 

O = C + vr p ( W - n") - V^VO, (4-4) 

in which n p and ix p as well as b pa are to be considered as independent free variables, while as 
before w pa is also to be considered as independent though not entirely free but given by (3.5) 



7 



which means that it is constrained to vary according to (3.9). Invariance of the ensuing action 
with respect to variation of the current n p itself gives the dynamical relation 

tt p = n p (4.5) 

from which the requirement of invariance with respect to the other allowed variations can be 
seen to give back the same dynamical equations as were given by the original Lagrangian L. 

Instead of considering the trivial relation (4.5) as a dynamical relation, let us now simply 
impose it as a defining restriction, so that the new Lagrangian (4.4) will reduce to the form 

£t = tt + &^(V^-^v) , ( 4 - 6 ) 

where 

^ = A-fi p n p . (4.7) 

It can be seen that the (fixed background) variation of this function \l/ will be simply expressible 

as 

5tf = -n p 8ix p - \\ pu bw pa , (4.8) 

which suggests that instead of regarding the components n p or equivalently N pup as independent 
field variables from which the components \i p or equivalently H pup are derived by partial differ- 
entiation according to (3.3), it would be convenient to treat the components fi p as independent 
variables, from which the components n p are obtained by partial differentiation according to 

(4.8) . 

This Legendre type transformation provides a new formulation in which the dual bivector 
b pa to the Kalb Ramond form comes in as a Lagrange multiplier imposing the condition that the 
vorticity 2- form should be the exterior derivative of the momentum 1-form. In this conjugate 
reformulation, \1/ is to be considered as a function just of the freely variable momentum covector 
fi p and of the independent but constrained vorticity 2-form w pa whose variation is governed by 

(3.9) . Under these conditions the corresponding variation of the new Lagrangian (4.6) will be 
given by 

6tf - V p {b p °5iia) = (W - n')5n P + (V p /i CT - l -w pa )5b p ° - l -(b p ° + \ p °)5w pa . (4.9) 

As well as giving back the relation (4.2) and the specification (2.3) of the vorticity as the 
exterior derivative of the momentum, the use of (3.9) in the new variational principle gives 
back the main dynamical relation (3.10) in the form 

fw pa = (4.10) 

where the modified current vector is given by the definition 

f = v ff (6^ + X pa ) = n p - V a \ ap , (4.11) 

from which it is apparent that, like the particle current vector n p itself, the modified current is 
also conserved: 

V p f = . (4.12) 
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Since the contribution to the action from the extra terms introduced in the transformation 
(4.4) does not depend on the metric, it follows that the stress momentum energy tensor derived 
from the modified Lagrangian C\ namely 

T p c = n p ii a - \ pv w va + , (4.13) 

will automatically be the same as the tensor (3.13) that was obtained from the original La- 
grangian C. 



5 Explicit development. 

The form of the master function A for the original formulation given in Section 3 can be made 
more explicit by expressing it in terms of the three independent scalars that can be constructed 
from the tensors n p and w pa on which it depends. These can be conveniently be taken to be 
the magnitude n of the particle current vector n p itself, the scalar magnitude w of the vorticity 
covector w pa and the magnitude ( of the associated Joukowsky lift force density vector ( p as 
defined by 

n 2 = -n p n p , w 2 = , ( 2 = ( P ( p , (5.1) 

where the Joukowsky vector is defined as 

Cp = w P «n° = l -wE p »N, vp . (5.2) 

This vector is interpretable as representing the volume density of force that would be exerted 
on the vortices as an expression of the Magnus effect, by the relative flux (if any) of the fluid 
according to the simple formula originally derived by Joukowsky (or Zhukovski, depending on 
how one transliterates from Cyrillic) for flow past a long aerofoil. The coefficients appearing in 
the derivation of A can be expressed explicitly in terms of partial derivatives with respect to 
these scalars in the form 

\H^ P = ~N^" - 3^wS^C p] , (5.3) 
r) A 8\ 

\ p ° = -2—w p ° - 4— . (5.4) 
The stress energy momentum tensor can thereby be rewritten in the manifestly symmetric form 

Tt = -^N^N™ - - 2^CC, + (A " 2^C 2 )< ■ (5.5) 

An analogously explicit analysis of the conjugate formulation developed in Section 4 can 
be made by expressing the generalised pressure function ^ in terms of the three independent 
scalars that can be constructed from the tensors \x p and w pa on which it depends. These can be 
taken to be the magnitude w of the vorticity vector, as already given by (3.8) or (5.1) together 
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with the magnitude p of the momentum covector and the magnitude h of the helicity vector as 
given by 

p 2 = -p p p p , h 2 = h p h p , (5.6) 
where the helicity vector' 16 ! is given by the definition 

h» = \e pup °p u w p(7 = wp p E™ , (5.7) 

which, by the degeneracy property (3.6), is such as to ensure that it will be conserved: 

V^/i" = . (5.8) 
In terms of partial derivatives with respect to w 2 , p 2 and h 2 one obtains the expressions 

no = 2—p p + 2 w pv>h a £«> , (5.9) 

\ p ° = -2—w pa - 2—h v H» p ° . (5.10) 
aw 2 oh 2 

The corresponding manifestly symmetric expression for the stress momentum energy density 
tensor (4.13) is found to have the form 

T ° = 2 d^ wPTw ™ + + 2 dV hPh * + ^ " 2 dh 2 ^ ' (5 ' U) 



6 Deviations from Perfect Fluid Limit. 



Not only in the usual terrestrial laboratory experiments but also in the neutron star applications 
for which the present theory is principally intended, the relevant macroscopic angular velocities 
are sufficiently small - and the ensuing microscopic vortex lines therefore sufficiently widely 
separated - for it to be a good approximation to consider the effect of the associated vorticity 
as a small perturbation to the dynamics. It will therefore be useful to formulate the theory in 
terms of deviations from a simple barytropic fluid model governed by an equation of state giving 
the zero vorticity density p just as a function of the particle number density n or equivalently, 
in conjugate, form giving the zero vorticity pressure P as a function of the effective mass per 
particle p. 

As a premilinary step, to prepare the way for a perturbation analysis in which the deviations 
will be considered to be small, we shall start without any loss of generality by considering 
deviations of unrestricted amplitude, simply decomposing the master function A as the sum of 
the perfect fluid contribution, p say, to which it reduces when the vorticity is zero, together 
with a vorticity dependent deviation term A in the form 

A = -p +A', P(d =p{ n }. (6.1) 
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We can also simply decompose the conjugate master function ^ in the analogous manner as 
the sum of the perfect fluid contribution, P e say, to which it reduces when the vorticity is zero, 
together with a vorticity dependent deviation term \1> in the form 

* = p e + tr , p e = pm . (6.2) 

With respect to the vorticity tensor w pa itself, the partial derivatives of these deviation terms 
will be given by the same tensor A pCT as for the total, but they will determine a reduced mo- 
mentum covector fi' p , and a reduced current vector n p according to the specifications 

SA' = ti' p 6n p -^X pff 6w f „ , (6.3) 

SV^-rCSnr-^Swp. . (6.4) 

Defining the particle current reference state unit vector u p ^ and the corresponding perfect fluid 
momentum contribution ji &p by 

n p = nu p , /x 0P = fi u 0P , // = /x{n} = ^ , (6.5) 
one obtains the total momentum contribution from (6.1) in the form 

Hp = H 0P + fip • (6.6) 

Similarly defining the momentum reference state unit vector u p and the corresponding perfect 
fluid particle current contribution n p Q by 

dP 

Hp = fiu eP , n p Q = n e u p Q , n e = n{fi} = , (6.7) 

one obtains the total particle current contribution from (6.2) in the form 

n p = n p + ri p . (6.8) 
Defining the particle current reference state perfect fluid contribution by 

T f = P.« + W + «l P , = - P. . (6-9) 
one can express the total stress momentum energy density contribution in the form 

T p a = T p a + T p a: (6.10) 

with 

T\ = n>; + \ pv w au + {K-n v n' v )gP a . (6.11) 

This can be interpreted as the stress momentum energy density contribution of a generalised 
kind of Stachel-Letelier type' 17 " 18 " 19 ' "string fluid", given by the vorticity dependent action 
contribution A'. 
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To obtain the analogous contribution from in the conjugate formulation we similarly 
define the momentum reference state perfect fluid contribution by 

T e = Pe« + P eGT + «)> Pe = We " P e ■ (6-12) 

We can then express the total stress momentum energy density contribution in the alternative 
form 

T^ = T^ + T\ , (6.13) 

with 

T\ = ri p fi a + y»w au + Vg" a . (6.14) 
It can be seen from (4.7) that the mutually conjugate deviation contributions will be related 

by 

A' - V = P e - P & + n v n'„ = p - p e + n>„ . (6.15) 

In the limit as deviations from the zero vorticity state tend to zero, it can be seen from the 
defining relations (2.1) that we shall have 

P e - P ~ n (P ~ Pe) > Pe ~Pe ~ P( n ~ n e) • ( 6 - 16 ) 

Since by (6.6) and (6.8) we shall have 

A* -A*~«aAV, n e - n ~ u eP ri p , (6.17) 

it can be seen from (6.15) that in this limit the deviations of the two mutually conjugate kinds 
of master function will coincide, i.e. we shall have 

^ S ~A\ (6.18) 

7 Asymptotically Separable Model for the Weak Vortic- 
ity Limit. 

A further simplification can be obtained by postulating that the master function has a sepa- 
rable form such that the vorticity dependence in the deviation term is contained in a negative 
factor, — T say, depending only on the scalar magnitude w, with a positive coefficient $ 2 that 
depends only on the relevant perfect fluid reference state: this state will be labelled by when 
determined by the particle current magnitude n so that the deviation A' defined in (6.1) will 
have the form A' = — $ 2 T; alternatively the reference state will be labelled by e when de- 
termined by the effective mass per particle \i so that the corresponding deviation defined by 
(6.2) will have the form = — $^T. It must be noticed that, although in general a model 
that is separable in the particle number representation (6.1) will not be exactly separable in 
the corresponding chemical potential representation (6.2) and conversely, it can be seen from 
the work of the preceeding section that the property of separability is nevertheless conserved 



12 



by the Legendre transformation in the weak vorticity limit. Such an asymptotically separable 
form, as given by 

A'~-$», ^~-$ 2 e TW, (7.1) 

in the relevant weak vorticity limit characterised by T — > 0, is suggested by our recent 
analysis^ 11 ' of the average stress momentum energy density for an individual vortex cell. The 
results of this analysis can be matched in a very satisfactory manner by an appropriate choice 
of the single variable functions T and $. 

Dropping the label or e in asymptotic formulae since there is no need to distinguish 
between the two kinds of reference state because either interpretation would be valid, as il- 
lustrated by the relation $ ~ $ ~ $ e , we can express the ensuing asymptotic forms of the 
particle current and momentum deviations given by (6.3) and (6.4) in the form 



and we similarly obtain 



H P ~ — — n P , n p ~ — n p , (7.2) 

[i an n an 

„ Xe^s x ~ $ 2 — . (7.3) 

2 aw 

On substituting this in the expression (6.11) for the deviation of the stress momentum en- 
ergy tensor from that of the particle current reference state, due to the vorticity dependent 
action contribution A', we see that the asymptotic form of this deviation will be expressible 
in terms of the fundamental tangential and normal projection tensors associated with the flux 
2-surfaces' 2 °]' 21 ], as given respectively by 



rf 9 = E p v E\ , (7.4) 



and 

in the form 



Af a = g? a - rf 9 , (7.5) 

in which, relative to the particle current reference state, the relevant effective vorticity pressure 
IT, acting transversely to the vortex flux lines, will be given by 

r/cb 2 

n<~ w x-(tf-n—)r . (7.7) 

In the equivalent conjugate formulation, the contribution from giving the deviation of the 
stress momentum energy density tensor with respect to that of the particle momentum reference 
state, will be analogously expressible in the form 

TZ ~ -/i— Tii'u, - $ 2 Ty; + iT L% , (7.8) 
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in which, relative to the particle momentum reference state, the effective vorticity pressure 77 s , 
relative to the momentum reference state, will be given by 



iT ~ w\ - $ 2 T . (7.9) 

In order to get the formulae (7.6), (7.7), (7.8) and (7.9) to agree with the corresponding 
formulae as obtained in the preceeding work' 15 ! by averaging over an individual vortex cell with 
circulation given by the standard quantum unit k = 2nh, it suffices to take the functions <3> and 
T to be given by 

$ 2 = ^, T = ^lnf^}, (7.10) 
fx 4 L w > 

where w Q is a fixed cut off vorticity value - whose exact value is unimportant in the limit 
w << w e under consideration. The specification (7.10) means that the function $ is to be 
identified with the dilation amplitude that has been shown to have a specially important role in 
the newly developed variational formulation of perfect fluid theory' 9 ! . The logarithmic derivative 
of this function is expressible in terms of the corresponding sound speed (2.2) according to either 
of the equivalent mutually conjugate specifications 

~ l - c s > 7^7777 = c s - 1 • ( 7 - n ) 



$ 2 dn s ' $ 2 dpi ~ s 

Since (7.10) also implies 

Aw~$ 2 T, (7.12) 

the deviation with respect to the particle current reference state is thus obtainable from (7.6) 
and (7.7) in the form 

T'? a ~ $ 2 T ((1 - c 2 )^ - C y a + (1 - c 2 ) lf a ) , (7.13) 

while the conjugate formula for the deviation with respect to the particle momentum reference 
state is similarly obtainable from (7.8) and (7.9) in the form 

TZ ~ $ 2 T((1 - a" V M(T - 7f CT ) , iT ~ . (7.14) 

The upshot of the foregoing analysis is that a master function given in terms of functions 
of the form (7.10) by 

A = -p -$ 2 TM +o{T}, (7.15) 

(where the label indicates that the quantities concerned are to be considered as functions just 
of the particle number density n) provides a vortex fibration model whose stress momentum 
energy density tensor is consistent with what has been derived' 15 ' by averaging that of an 
individual vortex cell in the weak vorticity limit. 

So long as w varies within a few orders of magnitude of some mean value (w) that is itself a 
great many orders of magnitude smaller than the the cut off value w Q appearing in (7.10), i.e. 

w w (w) , (w) << w e , (7.16) 
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which will be a very good approximation in typical contexts that can be envisaged for the 
application of the present theory, and in particular in the case of the neutron star matter for 
which it is principally intended, then the variation of the logarithmic factor in the formula 
(7.10) will be negligible, so that it will suffice to replace it simply by a linear function, with a 
constant coefficient /C proportional to the circulation k round an individual vortex, that will 
be given by 

T = !Cw, IC = -ln\^\ = hn\^\. (7.17) 
8tt l(w)J 4 l(w)J 



8 Dilatonic model. 



As an ansatz for matching the available results obtained' 15 ' from the analysis of individual 
vortex cells, the asymptotic form (7.15) obtained in the preceeding section is not the only 
possibility, but it is the the only one that satisfies the simplifying condition that the master 
function A be dependent only on the first two scalars, n and w, of the triplet introduced in 
(5.1), or equivalently that the conjugate master function \1/ depend only on w and the first 
scalar, /i, of the pair introduced in (5.6): 

<9A 

sT * aR = °- «»> 

Unlike an alternative simplifying assumption implicit in the earlier work of Lebedev and Kha- 
latnikov[8] which will be described in the appendix, the condition (8.1) has the convenient 
feature of treating the mutually conjugate functions A and \& on the same footing. The extent 
to which such a simplifying condition is accurate in the weak vorticity limit with which we are 
concerned remains to be checked by future work on non-axisymmetric vortex cells, but until 
there is any evidence that a more complicated ansatz may be needed, the most reasonable pro- 
cedure is to adopt what is obviously the simplest provisionally admissible supposition, namely 
(8.1), as a tentative working hypothesis. 

Proceeding on this basis, which means that the general laws (3.3) and (4.8) will simplify 
just to 

5 A = -ii 5n - A 5w (8.2) 

and 

5^ = n5fi- \5w (8.3) 

with 

= A + n/x , (8.4) 
the relations (5.3), (5.4) or their conjugate forms (5.9), (5.10), will reduce just to 

\p° = - w p° , (8.5) 
w 

and 

P = $-2 n P $-2 = » /g >6 N 

n 
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so that the expression (4.13) for the stress momentum energy density tensor will reduce simply 
to 

T p a = <$>~ 2 n p n a - \w' l w pv w va + ^g p a . (8.7) 

As a compromise between the alternative mutually conjugate master functions A and it 
is particularly convenient in such a case to work with a potential energy function U defined by 

W = -I(A + ¥), (8.8) 

whose variation law, subject to (8.1), is obtainable from (8.2) and (8.3) in the form 

SU = \5w+ ^5$. (8.9) 

This shows that - except in the degenerate "stiff" case for which fi is just proportional to n so 
that $ is constant and <5$ vanishes - the potential U can be used as a master function whose 
specification in terms of the independent variables w and <£> determines the original master 
function A and its conjugate \l/ by the formulae 

A = - — <5>- 2 - U , ^ = ^-$ 2 -W (8.10) 

where the quantities n 2 and ji 2 are obtained as functions of the vorticity amplitude w and the 
"dilatonic" amplitude $ via the partial differential relations 

9 ^ 1 9U n ^ndU . -i -i \ 

M = m ' n=$ V (8 - n) 

Subject to the assumption (8.1) the potential defined by (8.9) can be used to rewrite the 
Lagrangian (3.2) in the form 

C = - V V - V m X\ P X\ - w} , (8.12) 

which can be used as the basis of a new variational formulation in which the independent field 
variables are classifiable in three subsets: the first consists of the dual Kalb Ramond bivector 
components b pa , as defined by (4.1) or (4.2), in terms of which the current amplitude n is 
defined by 

n 2 = -(WK^ ; (8.13) 

the second subset consists of the vorticity base coordinates x 1 an d X 2 as defined by (3.5), in 
terms of which, with respect to a flat symplectic gauge (3.6), the vorticity amplitude will be 
given by 

w 2 = 2x) [ Xu ] (^X 1 )^X 2 ; (8.14) 

finally the third subset consists just of the dilatonic amplitude $ which is to be considered here 
as an independently variable scalar in its own right - except in the degenerate "stiff" case, for 
which U depends only on w with $ acting merely as a fixed coupling constant. 
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An attractive feature of the "dilatonic" formulation set up in this way is that it allows the 
physically required asymptotic form (7.15) to be matched by taking the new master function U 
to have a form that is not just separable but more specifically linearly separable with respect 
to <3> 2 , meaning that it is expressible as 

w } = V{®} + $ 2 T{u>} , (8.15) 

for suitable single variable functions V of $ and T of w. The condition of matching (7.15) in 
the weak vorticity limit leaves no freedom of choice in the specification of these single variable 
functions: T here must be the same as the function that is denoted by the same symbol in (7.15) 
and that is given explicitly by (7.10) or with sufficient precision for practical purposes by (7.17), 
while V must be the same as the function denoted by the same symbol in a recent discussion 
of the dilatonic formulation of the perfect fluid limit ^ , which means that it is determined by 
the equation of state for the mass density p of the underlying perfect fluid as a function of its 
conserved number density n according to the parametric specification 

V{n} = P {n}- 1 -n^, = n" 1 ^. (8.16) 

The exceptional Zel'dovichl 10 ] case characterised by a fixed dilatonic amplitude $ = a say arises 
from an equation of state of the "stiff" type p = n 2 /2a 2 + b where b is also fixed (acting just as 
a cosmological constant) which gives V = b and hence U = b + a 2 T. 

In any such linearly separable model the vorticity coefficient A will be given just by 

A = $ 2 ^ , (8.17) 
aw 

and, except in the degenerate "stiff" case (for which dV/d<& is indeterminate), the original 
master function will have a value that is given by a linearly separable expression of its own, 
namely 

A = _®dV _ V _ 2$ 2 T (81g) 

while the conjugate function \1/ turns out to have a corresponding expression in which the 
dependence on w drops out altogether, leaving just 

¥ = *^-V. (8.19) 

In applications for which the linear formula (7.17) for the vorticity dependent factor T is 
considered to be sufficiently accurate, so that (8.17) will take the specific form 

A = /C$ 2 , (8.20) 

while (8.11) will take the specific form 

n 2 = + , (8.21) 
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the Lagrangian (8.12) of the dilatonic formulation can be written out with a fully explicit 
presentation of the gradient dependence, which in this case will be homogeneously quadratic, 

as 

C = ^2 WW™ - V° X \ P X\ - /C<&V2 X ^X% ] (V^ 1 )V«'X a -V{<&} . (8.22) 

In the application of the corresponding action principle, variation with respect to the dual Kalb 
Ramond bivector components b pa gives the formula (2.3) for the vorticity form (3.6) in terms 
of the momentum covector given by (8.6) in conjunction with the definition (4.2); variation 
with respect to the vorticity base coordinates x 1 an d X 2 with the definition (8.5) gives the basic 
dynamical equation of motion (3.10); while finally - except in the "stiff" case for which $ is 
simply fixed as a coupling constant - variation with respect to the dilatonic amplitude itself 
gives the equation (8.21) whose solution (for a given form of the equation of state function 
V{$}) determines the value of this amplitude $ as a function of the current magnitude n given 
by (8.13). 

In simple cases the correspondence (8.16) between the equation of state for p as a function 
of n and the associated equation of state for V as a function of $ can be made explicit as 
discussed in the recent analysis of the perfect fluid case' 9 '. The simplest example of relevance 
as an approximation for neutron star matter is that of the standard relativistic polytrope with 
index 7 = 4/3 as given, for a constant coefficient a, by 

p = ^ n 4/3 ^ \/ = 2a 3 $ 4 . (8.23) 
2 

It is not to be expected that such a simple model could provide a completely accurate de- 
scription for the most general kinds of astrophysical and terrestrial applications. High accuracy 
would require an even more elaborate treatment allowing for the effect ' 22 1 of vortex lattice rigid- 
ity, and for many purposes in the context of neutron star applications it will also be necessary 
to allow for the further complication of magnetic effects involving the independently conducting 
proton fluid component W 2 ^. However the most serious limitation of the above model may turn 
out to be due to inadequate allowance for the effect of relative flow between the averaged particle 
current and the vorticity surfaces. Such an effect was deliberately excluded from consideration 
in the preceeding microscopic analysis' 15 ^ of an individual vortex cell that provided the averages 
to which the macroscopic model proposed here has been matched. A more difficult microscopic 
analysis (lacking the cylindrical symmetry that facilitated the previous calculation) still needs 
to be carried out by future work to determine how the relevant averages are affected by the 
inclusion of the effect of relative flow. Until the results of such a more complete microscopic 
analysis are available, the simple model set up in the present and preceeding section should 
be considered to be physically trustworthy only in the limit of negligibly small values of the 
Joukowsky vector (5.2) which can be considered as specifying the magnitude and direction of 
the relative flow. 

We may sum up by saying that within the subcategory characterised by the ansatz (8.1) 
postulating that the master function A depends only on the current magnitude n and the 
vorticity scalar w but not on the Joukowsky scalar (, the condition of matching the results of 
our preliminary microscopic analysis' 15 ' for the case £ = leads uniquely, as far as the weak 
vorticity limit is concerned, to a model of the kind set up in Section 7 and Section 8. A more 
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complete microscopic analysis may show that within the general framework provided in Section 
3 some other ansatz, such as the example® described in Appendix B, whereby A is made to 
depend also on (, will be necessary for matching cases for which ( differs significantly from 
zero. Nevertheless the obvious ansatz of ( independence has provided an elegant prototype 
model, which should provide a reasonably accurate description of stationary configurations in 
which there is no relative flow of the fluid relative to the vortex congruence, and perhaps even 
a moderately realistic description under more general circumstances. Considered as a toy, this 
prototype model can also be used, as described in Appendix A, to provide a more general 
illustration of the application of cosmic string theoretical concepts to superfluid dynamics that 
was considered by Davis and Shellard' 4 '® within the restrictive framework of the "stiff" limit 
characterised by a fixed value of $ with vanishing V. 
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Appendix A: The tie up with cosmic string theory. 

The preceeding theory is designed to represent the macroscopically averaged effect of a 
congruence of microscopic vortices of which each individual member can be considered as an 
example of the kind of topological defect known in a cosmo logical context as a global cosmic 
string, of which the simplest type is provided by the axion field theory in which the relevant 
underlying fluid model is of the "stiff" kind in which the dilatonic amplitude <3> has a constant 
value. (The qualification "global" is needed here to distinguish such a relatively extended 
configuration from the more strictly string-like gauge coupled case for which, instead of being 
logarithmically divergent, the defect energy distribution is locally confined: such "local" string 
defects occur in the familiar laboratory context of ordinary metallic superconductors, and are 
also predicted - in coexistence with the dynamically dominant global string defects with which 
the present work is concerned - in the interior of neutron stars' 1 '®, where the relevant gauge 
coupled superconducting current is constituted from protons.) 

To see how the theory set up in Section 8 is to be interpreted in the terms of the technical 
machinery originally developed to describe string defects of the cosmic type, it is to be remarked 
that the general form (8.7) that is obtained from (8.1) for the stress momentum energy density 
tensor is decomposable in a natural manner as the sum, 

n = 7* - ^Trf , (Al) 

of a spacially isotropic (non-barytropic) perfect fluid type contribution 

T p a = n"n a + (* + Xw)gf a (A2) 

together with a term of the form that would be given by a Stachel-Letelier' 17 ^ 18 ^ 1 ® type "string 
fluid", whereby each individual vortex cell contributes as if it were a string of the simple 
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degenerate (longitudinally Lorentz invariant) Goto-Nambu type with fundamental tensor r\ p a 
and tension T given by 

T = k\, (A3) 

since the vorticity w is interpretable as due to a flux of w/k distinct microscopic vortices per 
unit area, where k denotes the constant momentum circulation associated with each individual 
vortex. 

More particularly, for the specific model of the linearly separable kind characterised by 
(8.15) with the vorticity dependent factor T given by the linear formula (7.17) on which (8.22) 
is based, the tension of each such (degenerate Goto-Nambu type) string will be given by 

T = /t O> 2 = — $ 2 ln|^} . (A4) 

and in this case the perfect fluid contribution (A2) will be expressible simply by 

T£ = + (\n^ - V)g\ (A5) 

which has just the form that would be obtained from (8.22) if the vorticity coupling coefficient 
K, in it were set to zero. 

Since the closure condition (3.7) gives the kinematic identity 

V u (w£ vp ) =0 w pa V u w V(T = w A. v p V u w - w 2 K p , (AQ) 

where K p is the geometric curvature vector of the vortex sheet, as given' 21 ] in terms of its 
fundamental tensor (7.6) by 

K p = r? c V v rf P = , (A7) 

the proportionality relation (8.5) can be used to reduce the dynamical equation of motion (3.10) 
to the form 

C p = w(V p V v \-\K p ) (AS) 

for any model of the subcategory characterised by (8.1). Since the vorticity w is interpretable 
as due to a flux of w/k distinct string-like vortices per unit area, the Joukowski force density 
( p , as defined by (5.2), can be interpreted as representing a corresponding average Joukowski 
force per unit length, Z p say, acting on each individual vortex cell, that will be given by 

Z P = -C P = \^ u N, vp . (A9) 

The equation of motion (A8) can thus be seen to be expressible as the condition that this 
Joukowski force per unit length should be given in terms of the corresponding effective vortex 
tension (A3) simply by 

Z p + TK p =r p V p T . (A10) 

This can be recognised as a generalisation of the standard equation of motion for a global cosmic 
string in an axion field background' 13 ^ 14 ], to which it reduces when the tension T is constant so 
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that the gradient term on the right drops out. The constant tension case arises for the model 
(8.22) obtained from (7.17) to which (A4) applies when the equation of state is of the "stiff" 
type for which $ is constant. 



Appendix B: The Lebedev Khalatnikov subcategory. 

It is easy to confirm that the category of theoretical models obtained by the convective vari- 
ational approach described in Section 3 includes as a subcategory the class of models originally 
set up by Lebedev and Khalatnikov' 8 ' . To show this it suffices to identify their generalised pres- 
sure function \1/ with the Legendre transformed master function ^ on which the reformulated 
variational principle of the Section 4 is based, since on this basis our equation of motion (4.10) 
agrees precisely with the one that they obtained using a Clebsch type variational procedure. 
Instead of using a Lagrangian of the form (4.6) in which the master function is supplemented 
by the appropriate Kalb Ramond type coupling term, the variational formulation proposed by 
Lebedev and Khalatnikov used ^ by itself as Lagrangian, the cost of this apparent simplifica- 
tion being the need to consider the relevant momentum covector \i p not as an independent field 
variable in its own right but merely as a derived quantity, specified by an appropriate set of 
Clebsch type potentials and their gradients. The question of whether one prefers to work with a 
set of scalar Clebsch type potentials as in the approach pioneered by Lebedev and Khalatnikov, 
or with a single antisymmetric tensorial potential of Kalb Ramond type as in the approach 
used here is merely a matter of taste and convenience in view of the ultimate agreement of the 
ensuing field equations. The Clebsch formulation has the technical advantage of using fewer 
independent field components but has the drawback that a larger number of different alphabet- 
ical letters are needed to characterise them. From the point of view of mathematical elegance 
and ease of formal manipulation, the Kalb Ramond formulation (either in its original version 
as given in the section 3 or the dual version of section 4) is more satisfactory: it involves a 
greater degree of gauge dependence but has the advantage of being more economical in its use 
of algebraic symbols. 

Although the postulates of the Lebedev Khalatnikov theory are thus confirmed to be con- 
sistent with those adopted here, their category of models was restricted to a subclass of those 
considered in Section 4 by the supposition that the master function \1> depended only on ji and 
h but not on the third independent scalar w, i.e. they assumed that the form of the master 
function would be characterised by 

Subject to this restriction the variation of \I/ would be given simply by 

5^ = £ p 5h p - m p 5fi p , (B2) 

with 

i f = 2-k f , = (B3) 

In terms of these quantities, the stress momentum energy density tensor (5.11) will be reducible 
to the simple form 

T% = w?n a + hn a + (* - l v h v )g\ , (54) 
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which agrees with the expression given by Lebedev and Khalatnikov' 8 ' . However a small dis- 
crepancy does occur in their expression for the modified current vector (4.11), which omitted 
the final term of the formula 



j p = m 1 



p _ H pav VJ v + 2wi v £ vp 



(B5) 



that we obtain in this case. 

The choice of the particular kind of model introduced in Section 7 and Section 8 to match 
the results of our recent microscopic analysis^ 15 ' for the case ( = was determined by the 
obvious ansatz that the master function A should not depend on the Joukowsky scalar (, but 
only on n and w. We saw that in the weak vorticity limit under consideration, this ansatz 
was equivalent to the postulate that the conjugate master function \1/ should depend only on 
H and w but not on the helicity scalar h. However there is not any obstacle to matching the 
same results^ 15 ' by a model characterised by the alternative Lebedev - Khalatnikov ansatz to 
the effect that \I> should instead depend on \x and h but not on w. This can be done simply by 
replacing w in the deviation term given by (7.1) by the ratio h//i which is the same as w when 
( vanishes. It can be checked directly that the model specified by the deviation term 



that is obtained in this way will have a current vector and a stress momentum energy density 
tensor that agree with those obtained from the ansatz adopted in Section 7 for small w provided 
( vanishes. However this ansatz leads to a more complicated form for the conjugate formulation 
which can be seen to have an explicit ( dependence given by 



Until further evidence is provided by a more complete microscopic analysis, the criterion of 
simplicity would seem to rule against this latter model in favour of the separable model proposed 
in Section 7. However it may well turn out that neither alternative is adequate for cases in 
which ( is too large to be neglected. 
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